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inclusive DIS at high energies. We present new results using the 1994 data 
from HI, and confirm directly the expected running of a^. 



Determining the strong coupling as from fixed target inclusive deep in- 
elastic scattering^ requires very precise measurements of structure functions 
at large x over a wide range of in order to extract the relatively small 
violations of Bjorken scaling. In colliding beam deep inelastic scattering ex- 
periments, such as are now being performed at HERA, structure functions may 
be measured at small x, where Bjorken scaling violations are much stronger. 
Indeed the steep rise in ^2(0;, Q^) as l/x and both increase is driven es- 
sentially by the thrcc-gluon vertex'^, resulting asymptotically in a new double 
scaling behaviour"^. The predicted slope of the rise in F2 has been confirmed 
very precisely using the recent HERA data*~^, and two loop corrections are 
now discernible Since the three gluon vertex is itself directly proportional 
to as, the success of double scaling immediately suggests^ that F2{x, Q"^) must 
be rather sensitive to the value of as in the double scaling region, and indeed 
this expectation is confirmed by simple double scaling fits to the data which 
include as as a free parameter^^'^. 

Of course a proper determination^^ of as from F2 at small x requires 
a full two loop calculation which matches on to all the other data (and in 
particular structure function and prompt photon data) at large x. Care must 
also be taken to disentangle the effects of higher twists and higher logarithms 
(of 1/x): fortunately this turns out to be rather easier than at large x, where 
the effects of both higher twists and infrared logarithms (of 1—x) may become 
very significant. Here we will briefiy review how such a small x determination 
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Figure 1: The three fitted parameters (A,, Xg,as{mz)) for a range of values of Qo- 
The distribution is fitted and evolved in MS. 




Figure 2: Contour plots of in the three orthogonal planes (A,, Ag), (as{mz), A,) 
and (as(»nz), Ag) through the global minimum. The first eleven contours are at 
intervals of one unit, while those thereafter are at intervals of five units. The 
starting distribution is in MS at Qo = 2GeV. To compute the statistical and 
systematic errors are combined in quadrature, and the normalization uncertainties 
included: at the minimum = 80 for 166 degrees of freedom, the normalizations 
are within 1% of the experimental values. 



may be performed, and then present some preliminary results obtained using 
the most recent data^ from the HI collaboration. 

Our basic procedure^^'^^ is to take a set of globally fitted parton distri- 
butions, for example those of CTEQ-'^^ or MRS^^, evolve them (using a two 
loop evolution code) to some new starting scale Qq, and there cut off the tails 
of the sea and gluon distributions, replacing them with new tails xq ~ x^" , 
xg ^ X for X ^ 0.01.* These new distributions, together with the valence 
distributions, are evolved up (or down) to the HERA data, again at two loops, 
using a particular value of as(m^). The three parameters (Ag, ^g^asi^z)) are 
then adjusted to minimise the to the data. The final distributions are by 

* lu practice this is achieved by refitting the large x distributions at Qo with \ and Ag 
kept fixed. 
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Figure 3: The strong coupling measured independently at five different scales. The 
horizontal bars indicate the width of each bin. The curve shows the theoretical 
two loop running when as(Tnz) = 0.122 (the sixth point). The points are over- 
correlated due to correlated systematics. 

construction a reasonable fit to all the data with x ^ 0.01 used to determine 
the original input distributions, since QCD evolution is causal in x. However 
the shape of the distributions in the HERA region is largely independent of the 
details of the input, since asymptotically it takes a universal (double scaling) 
form which depends rather sensitively on the value of a^.^ Higher order effects 
due to choosing large x distributions which have been fitted using different 
values of should in principle also be included. In this sense our procedure 
is similar to that used to determine ag from, for example, the 2 + 1 jet rate at 
HERA, or the inclusive jet rate at the Tcvatron. 

The results of such fits to the HI 1994 data* for a range of values of Qo are 
presented in fig. 1. The stability of the fitted value of ag in the region where 
the is lowest (2 GeV ^ Qq ^ 4 GeV) is a useful test of the effectiveness 
of the parametrization employed at Qo- More information about the shape of 
the minimum in the at a particular value of Qo is displayed in the contour 
plots fig. 2. The third plot shows that although the shape of the inferred 
gluon distribution and as are indeed correlated, this does not lead to a large 
uncertainty in as because of the close relation between the quark and the 
gluon shown in the first plot: very steeply rising or falling gluon distributions 
at 2 GeV are incompatible with the double scaling seen in the HI data. This 
explains why it is possible to determine ag at small x even though the gluon 
distribution is large there. 

Our preliminary result using the HI 1994 data is as{mz) = 0.122 ± 
0.004(exp). The theoretical error has not yet been determined, but should 
be a little less than that found^^ in our analysis of the 1993 data, partly be- 
cause the data themselves now severely limit the allowed size of corrections due 
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Figure 4: Contributions to the total from various structure function data sets 
included in global fits at six different values of as{Mz), adapted from ref.l4. The 
data are from BCDMS F2 (open squares, 142 points), NMC F2 (open diamonds, 74 
points), CCFR F2 (open circles, 80 points), CCFR F3 (open triangles, 80 points), 
and HERA F2 (asterisks, 93 points from HI + 56 from ZEUS, both from the 1993 
run). The curves are also from fits to HERA F2 data: the dotted (A^, = \g) and 
clashed (A^ ^ Ag) curves using 1993 data (122 points) are from the analysis in 
ref.l3, while the solid curve is from this analysis (HI 1994 data, 169 points). In all 
these curves as{Mz) in the large x input distribution was kept fixed. 

to higher order logarithms of l/.x.^" The error due to higher twist corrections 
seems also to be very small, since their size is again limited by the data.^° Fur- 
thermore the allowed range of variation of renormalization and factorization 
scales can now be limited by requiring that the fit to the data is not unduly 
worsened. A more complete analysis will be made when all the ZEUS 1994 
data have been published. 

By only including some of the data in the fit, over a limited range of Q^, we 
can actually perform independent measurements of as at different scales. The 
results of such a procedure are displayed in fig. 3. It is gratifying to observe 
directly the running of as in a single experiment. Of course we already knew 
that this had to work: the slope of the rise of F2 is in itself a direct measurement 
of the first coefficient of the /j-function.^ 

It is important to emphasise that our determination uses only HERA data 
directly: it is not a global fit. Some of the advantages and disadvantages of the 
two techniques may be seen by comparing our results to those of the subsequent 
MRS global fits described in ref.l4 (see fig. 4), or indeed to very similar results 
from CTEQ.^^ In a global fit the for each data set are in general simply 
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added together, irrespective of the fact that different experimental collabora- 
tions may treat their errors rather differently. The data from BCDMS then 
dominate the determination of as not because there are many more points, but 
because the per point turns out to be relatively large (a better fit to these 
data, with more parameters, was obtained in ref.l). Indeed the minima in the 
of the NMC and CCFR F2 data actually disappear in the MRS global fit. 
The minimum in the 1993 HERA data is uncorrupted, however, and is indeed 
consistent both in position and width with the direct analysis of ref.13: this 
is to be expected as the shape of F2 at small x is largely independent of the 
detailed structure at large x.^ Note that the global analysis^^ has — Xg-. 
relaxing this constraint softens the minimum but does not eliminate it. The 
present analysis of the 1994 data shows that in the future HERA structure 
function data will begin to play an increasingly significant role in the global 
fitting of as (although our solid curve in fig. 4 will presumably not fall quite 
so steeply once as is varied in the input distribution) . It is interesting that the 
relatively high value of as favoured by HERA is also found in analyses of the 
CDF and DO inclusive jet data.i^^i^.iv 

In conclusion, F2 at small x and large is an excellent place to measure 
as', the data is now very precise, the dependence on as is strong, one need only 
do simple fits with a small number of parameters, while on the theoretical side 
uncertainties due to higher logarithms seem to be unexpectedly small, while 
higher twists are truly negligible. Furthermore one can perform a direct test 
of the running of as at spacelike in one experiment. We hope that soon 
experimentalists will perform their own determinations, taking properly into 
account the correlations of their systematic errors. 
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